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ABSTRACT 

In order to characterize Damped Lyman Alpha systems (DLAs) potentially host- 
ing first stars, we present a novel approach to investigate DLAs in the context of 
Milky Way (MW) formation, along with their connection with the most metal-poor 
stars and local dwarf galaxies. The merger tree method previously developed is ex- 
tended to include inhomogeneous reionization and metal mixing, and it is validated 
by matching both the Metallicity Distribution Function of Galactic halo stars and 
the Fe-Luminosity relation of dSph galaxies. The model explains the observed Ahi-Fc 
relation of DLAs along with the chemical abundances of [Fe/H]< —2 systems. In this 
picture, the recently discovered z a b s as 2.34 C-enhanced DLA (Cooke et al. 2011a), 
pertains to a new class of absorbers hosting first stars along with second-generation 
long-living low- mass stars. These "PopIII DLAs" arc the descendants of RVcooling 
minihalos with Mh ~ 10 7 M Q , that virialize at z > 8 in neutral, primordial regions of 
the MW environment and passively evolve after a short initial period of star forma- 
tion. The gas in these systems is warm T g ps (40 — 1000) K, and strongly C-enriched 
by long- living, extremely metal-poor stars of total mass M* as 10 2-4 M Q . 
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1 MOTIVATION 

Damped Lya absorption systems (DLAs) are high-column 
density neutral gas reservoirs, Am > 10 20 ' 3 cm -2 , observed 
at intermediate redshifts, z < 5, in the spectra of dis- 
. tant quasars. Although their nature is still unclear (Pet- 
tini 2004), the key role of DLAs to understand galaxy 
formation (Wolfe, Gawiser & Prochaska 2005) is widely 
recognized. So far, more than 1000 DLAs have been ob- 
served and the iron abundance measured in w 150 sys- 
tems [Fe/H] « [-3.5, -0.5] (Prochaska et al. 2007). Among 
these, very metal-poor (VMP) DLAs with [Fe/H]< -2, can 
be used to study the initial phases of heavy element en- 
richment of the interstellar medium (ISM) of early galaxies. 
Indeed, if VMP stars observed today in the Galactic halo 
and in nearby dwarf spheroidal galaxies (dSphs) are the liv- 
ing fossils of the first stellar generations, VMP DLAs may 
well constitute the gas reservoir out of which such pristine 
stellar populations formed. Following the medium resolu- 
tion study of VMP DLAs by Penprase et al. (2010), Cooke 
et al. (2011b) have recently presented a high spectral res- 
olution sample, including 22 VMP systems. In these DLAs 
[C/0]« -0.3 and [Fe/0]w -0.4 independently of [Fe/H] 
and with little scatter, in agreement with measurements in 



VMP Galactic halo stars (Fabbian et al. 2009). So far, the 
only exception to this general trend is represented by a DLA 
with [Fe/H]~ -3 and N m = I0 2a55±0 10 cm" 2 observed in 
the spectrum of the QSO J0035-0918 at z aba « 2.34 (Cooke 
et al. 2011a). This system has [C/Fe]= 1.53, i.e. « 20 times 
larger than any other DLA. Moreover, its abundance pat- 
tern shows a clear 'odd-even' effect and is consistent with 
the predictions for the yields of Z — faint supernovae (SN) 
with m„ ~ 25A/q (Kobayashi et al. 2011). Are we observing 
for the first time a DLA whose gas retains the imprint left 
by the first stars? 

To characterize DLAs potentially hosting the first stars 
or their ashes, we propose a novel approach that simultane- 
ously follows the evolution and chemical properties of DLAs 
and their connection with the most metal-poor stars and 
galaxies observed in the Local Universe based on the re- 
sults obtained using the merger-tree code GAMETE (Sal- 
vadori, Schneider & Ferrara 2007, SSF07; Salvadori, Ferrara 
& Schneider 2008; Salvadori, & Ferrara 2009, SF09). 



2 MODEL SUMMARY AND VALIDATION 

The model key points can be summarized as follows (see 
SSF07 and SF09 for details): 
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(i) Hierarchical merger histories of a MW-sized dark mat- 
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Figure 1. Left panel: observed (points with error-bars) and simulated (histogram) metallicity distribution function of Galactic halo stars. 
The data points are the sample by Beers & Christlieb (2006) with the inclusion of the three hyper-iron poor stars (Christlieb et al. 2002, 
2006; Frebel et al. 2005). The histogram is the average MDF value over 50 merger history of the MW, re- normalized to the total number 
of observed stars with [Fe/H]< 2. The shaded area represents ±1ct errors. Right panel: observed (points with error-bars) and simulated 
(contours) iron-luminosity relation for the MW dSph galaxies. The data points are by Kirby et al. (2008) with the inclusion of the two 
faintest dSphs (Willman et al. 2005, Geha et al. 2009). The colored shaded areas correspond to regions that include the (99,95,68)% of 
the total number of possible dSph candidates selected in 50 merger histories of the MW. 



ter (DM) halo are reconstructed from z — 20 via a Monte 
Carlo algorithm based on the Extended Press-Schechter the- 
ory (SSF07). 

(ii) Star formation (SF) is followed along the tree in halos 
exceeding a mass threshold, M s f, whose evolution (Fig. 2) 
is governed by: (a) the photo-dissociating Lyman- Werner 
(LW) background, quenching Hb-formation in T V i r < f0 4 K 
minihalos; (b) the gas temperature in ionized regions of the 
Galactic Medium (GM), preventing gas-infall in halos with 
T V i r lower than a threshold value, T t h- We assume that SF 
is active in T„j r > T±h — 2 x 10 3 K minihalos at z > fO 
(Dijkstra et al. 2004). At lower z, T t h (and hence M s f) is 
assumed to linearly increase up to the value set by the end 
of reionization T th R! 2 x 10 4 K (Kitayama et al. 2000) for 

Z *C Zyei — 6. 

(iii) Inhomogeneous reionization is modeled by random 
sampling the reionization history implied by M s f(z) to 
switch off (on) gas accretion in minihalos that form in ion- 
ized (neutral) regions. 

(iv) The SF rate is taken to be proportional to the mass 
of cold gas, M = e*M g /tff, where e* is the SF efficiency 
and tff the halo free-fall time. In minihalos e„ is reduced 
as e = e»[l + (T vir /2 x 10 4 K)~ 3 ] -1 due to the ineffective 
cooling by H2 molecules (SF09). Low- mass, Population II 
(PopII) stars form according to a Larson IMF when the gas 
metallicity exceeds Z CI = 10~ 3 S Z Q (Schneider et al. 2002). 
At lower metallicity, PopIII stars form with a reference mass 
value m* = 25M© and explosion energy Esn = 10 51 erg 
consistent with faint SNe. 

(v) The abundance evolution of different chemical 

elements^ (from C to Zn) is traced in both the ISM and 

t For m* < 8Mq stars we use yields by van den Hoek & 
Groenewegen (1997) (Z > 10~ 3 ) and by Meynet & Maeder 
(2002) (Z < 10~ 5 ); for more massive stars we use Woosley & 
Weaver (1995) with a systematic halving of the Fe yield (Timmes, 



in the GM by taking into account mass- and metallicity- 
dependent stellar evolutionary timescales (Raiteri, Villata 
& Navarro 1996) and SN feedback (Salvadori, Ferrara & 
Schneider 2008). 

(vi) To account for the incomplete mixing of SN ejecta 
within the ISM of gas-poor galaxies, gas outflows have a 
metallicity Z m = Z ISM + r)(M g )Mz / 'M ej , where M z is the 
mass of newly formed metals, M e j is the mass of gas ejected 
out of the halo, and r\ is a function of the gas mass r\ = 
0.5+0.65 tanh[(loglO(M ff )-7.0)/2.0] that varies in the range 
r, = [0,1]. 

(vii) The probability for newly formed halos to reside in 
a GM metal enriched region is P(z) = Qz /Qs>s c , where 
Qz(z) = 1 — exp(T,i4TvRl(i)/VMw(z)) is the filling factor 
of metal bubbles within the MW physical volume^, and 
Qs>s c (z) is the volume filling factor of fluctuations with 
overdensity above the critical threshold, 5 > S c — 1.686, 
for the linear collapse (Miralda-Escude, Haehnelt & Rees 
1999). The latter quantity describes the abundance of high- 
density regions, in which metals first penetrate (Tornatore, 
Ferrara & Schneider 2007). Objects in enriched (primordial) 
regions are assigned an initial metallicity Z V i r — Zqai/Qz 
(Z v i r = 0). 

The model is calibrated by best-fitting the SF and feed- 
back efficiencies to reproduce the global properties of the 
MW (stellar/gas mass and metallicity). In Fig. 1 we show 
that the average metallicity distribution of [Fe/H]< —2 stars 
over 50 possible MW merger histories matches the observed 
Galactic halo MDF. In the same Figure, the theoretical Fe- 
Luminosity relation for dSph galaxies has been obtained by 
evolving in isolation star- forming halos with Mh < M2 a , 
candidates to remain MW satellites (Diemand, Madau & 

Woosley & Weaver 1995); yields for faint SNe are from Kobayashi 
et al. (2011). 

$ We assume V M w(z) = 5(1 + z) -3 Mpc 3 . 
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Figure 2. Dark matter content of MW satellites (Mh < M2a) as 
a function of their formation redshift. The two lines show the evo- 
lution of: the minimum mass of star-forming haloes, M s f (solid), 
the halo mass corresponding to 2a density peaks (dashed). The 
dotted shaded area identifies the region populated by H2-cooling 
minihalos with T v i r < 10 4 K. The number labels and colors 
specify different populations of satellites (see the text). The gray 
shaded area delimits the evolution of M 3 f predicted by assuming 
an early (upper limit) and a late (lower limit) reionization history 
(Gallerani et al. 2007). 



Moore 2001; SF09 for details). In Fig. 2 M sf (z) is compared 
with that obtained by using data-constrained reionization 
models. To this end we compute the dissociating LW back- 
ground intensity associated with an early /late reionization 
history (Gallerani, Choudhury & Ferrara 2006) by using the 
results provided in Fig. 6 by Ahn et al. (2009); we then con- 
vert the flux into a critical mass following Machacek, Brian 
& Abel (2001) and extrapolate to higher masses. At each 
redshift our M s f is within the range allowed by the two 
reionization histories (gray shaded area). 



3 THE DWARF GALAXY ZOO 

We now apply our model to DLAs starting from the ori- 
gin of the C-enhanced, [Fe/H]« -3 DLA at z abs = 2.34. 
Two possibilities must be considered: (a) a recently viri- 
alized halo, later merging into the MW, or (b) a satellite 
galaxy which assembled at earlier epochs and evolved in 
isolation. The first possibility can be excluded on the ba- 
sis of a metallicity argument: metal- free regions of the MW 
environment are expected to disappear at (z) ~ 8.7 ± 1.5 
when P — Qz /Qsyi.asa > 1- Afterward, newly vitaliz- 
ing halos form from metal enriched GM gas, whose initial 
abundances [X/H]„«, = ([X/H])gm /Qz ■ Since at z — 2.3, 
Qz = 0.989 ± 0.005 and ([Fe/H]) = -1.28 ± 0.02, proto- 
galaxies with an initial iron abundance [Fe/H]< —1.3 cannot 
form anymore. 

In the second hypothesis, the dark matter (DM) halo 
hosting the DLA must be associated to density fluctuations 
< 2a at its final assembling epoch in order to evolve in iso- 
lation (no further merger or accretion) and become a satel- 
lite. By inspecting Fig. 2 we can distinguish among different 
populations of objects that satisfy this condition: (1) Lya- 
cooling halos that assembled after the end of reionization; 



Fornax and Leol, the most luminous among the observed 
dSphs, are predicted to belong to this domain; (2) star- 
forming Lya-cooling halos that assembled before the end of 
reionization; classical dSphs, L > 1O 5 L0, such as Sculptor 
(SF09) are members of this population. (3) Eb-cooling, inef- 
ficiently star-forming minihalos that appeared at z ~ 8 — 10; 
this is the domain of ultra-faint dSphs (SF09). (4) "sterile" 
halos that formed before the end of reionization and unable 
to trigger SF (Mh < M s f). We selected the satellite candi- 
dates in 50 possible merger histories of the MW, and evolve 
them in isolation down to redshift z = 2.34, when we com- 
pute their H i column density, Am, assuming that the gas 
is fully neutral with a molecular weight n = 1.22: 
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the gas radius is written as r g = ar v i r , with a = 0.18 
for star-forming galaxies, that have presumably developed 
a disk (Ferrara, Pettini & Shchekinov 2000), and a = 0.8 
otherwise. DLAs are canonically defined as systems with 
N m > 2 x 10 20 cm" 2 (Wolfe et al. 2005). Interestingly, we 
find that systems which populate Zone (2), i.e. the progeni- 
tors of classical dSphs, have already exhausted most of their 
gas by z = 2.34, i.e. log(A?m/cni~ 2 ) < 20.3. On the other 
hand, many objects among those formed in Zones (1), (3) 
and (4), have H i column densities compatible with DLAs. 
We will then focus on these candidates. 



4 DLAS AND FIRST STARS 

A comparison between the properties of candidate and ob- 
served DLAs at z ~ 2.34 is given in Fig. 3 and Fig. 4. Ob- 
jects formed in different Zones result in distinct DLA popu- 
lations, differing in their chemical abundances and Njjx- Sys- 
tems in Zone (3) are segregated in the C-rich, metal-poor, 
low- Am regions of the plots, matching the properties of the 
C-enhanced DLA. What is the origin of these systems? 

According to our model, C-enhanced DLAs are hosted 
by Mh W 6 ' 9 ~ 7 ' 2 Mq H2-cooling minihalos, which vital- 
ized before the end of reionization out of metal-free, neu- 
tral regions of the MW environment. Because of their low 
DM content the star formation history of such a primordial 
proto-galaxies is extremely short. Some tens of Myrs after 
the onset of star-formation the intensity of LW background 
becomes high enough (M s f > Mh) to suppress further SF 
activity, thus turning them into "sterile" minihalos. Metal- 
free stars are hosted by this population of DLAs, to whom 
we will hereafter refer as PopIII DLAs. 

The metal-enrichment by first stars gently proceeds in 
these objects, mainly because of ineffective cooling by H2 
molecules and their low gas-mass content, M g 
,6,3 Mq, favoring metals and gas loss (M^ « 



10 b 



io 5 - 6 m q ). 

As soon as Z > Z cr = 10~ 3 8 Zq, however, "normal" PopII 
stars can form and contribute to enrichment. Given the 
yields by faint SNe, the PopIII-to-PopII transition occurs 
when [Fe/H] cr ~ —4.8, implying that most of the iron ob- 
served in these systems originates from these second gener- 
ations of stars. The mass of relic stars in PopIII DLAs is 
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Figure 3. Gas abundances of DLAs with respect to the solar values (Asplund et al. 2009). Points are the observed values in VMP DLAs 
at z ab3 = (2.07 — 4.21) (triangles, Cooke+2011b) and in the C-enhanced DLA at z abs = 2.34 (square, Cooke+2011a). The color shaded 
areas show the results of the model for the three different populations of DLAs at 2 fx 2.34: PopIII (blue), very metal-poor (violet), and 
metal-poor DLAs (orange). For each population the intensity of the colors correspond to regions containing the (99, 95, 68)% of DLAs. 



expected to vary between M* m 10 2 ~ 4 Mq, the most star- 
rich systems being the most enriched ones. The rise of the 
C abundance at increasing metallicity (Fig. 3) reflects the 
gradual contribution by low-metallicity AGB stars, mainly 
producing C, some O, and limited amounts of N (for sys- 
tems with [Fe/H]= -3.0 ± 0.2 we find [C/H]= -1.6 ± 0.5, 
[N/H]= -3.8 ±0.9 and [0/H]= -2.4 ±0.4). The associated 
A^hi decrease (Fig. 4) is caused by gas consumption in the 
most star-rich systems, due to both astration and gas loss. 
Note that our model does not predict any [C/O] increase to- 
wards low [O/H] values, as reported by medium resolution 
observations (Penprase et al. 2010). 

The properties of VMP DLAs are well reproduced (Fig. 
3) by [Fe/H]< —2 starless systems formed in Zone (4), pas- 
sively evolving since their assembly epoch, 6 < z < 10. They 
form through merging of primordial and metal-enriched pro- 
genitors, virialized out of neutral regions of the MW envi- 
ronment. The metallicity spread of these systems depends 
on the gas enrichment in progenitor minihalos. Abundance 
ratios, instead, closely reflect those of the GM at the time 
of formation. Although the overall GM metallicity increases 
with time, abundances ratios get locked to the dominant 
stellar population which contribute to the enrichment, i.e. 
type II SNe. As a consequence, [C/O] and [C/Fe] ratios 
show little dispersion and are independent of Z, but also 
of redshifts, because of the passive evolution of these star- 
less DLAs. These findings are in perfect agreement with the 
new observational results by Becker et al. (2011). 

Finally, [Fe/H]> -2 DLAs are T vir > 10 4 K halos form- 
ing through merging of star-rich progenitors either before 
(Zone 4) or after (Zone 1) the end of reionization. Their sys- 
tematically higher Z is a consequence of self-enrichment by a 
substantial (M* ~ 10 6 — 10 9,5 Mq) population of long-living 
stars. Similarly, the high [C/Fe] ratios reflect the strong con- 
tribution by AGB stars. The gas mass in these DLAs is 
M g = 10 7 - 1O 9 - 5 M , thus resulting in higher N m values, 
despite of the lower formation z and larger Mh ~ 10 8_11 Mq 
(eq. 1). Note however that the A^hi derived for these DLAs 
must be interpreted as an upper limit. In fact, Mh > M B f 
objects (Fig. 3) are actively star-forming at z = 2.34, with 
M fa (0.1 — lO)M0yr -1 , and hence part of their gas is pre- 
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Figure 4. Ajji distribution vs. iron abundances of DLAs. Points 
with error-bars are observed values for metal-poor DLAs at 
tabs = (1-9 - 4.9) (circles, Prochaska et al. 2007), VMP DLAs 
at z abs = (2.07 - 4.21) (triangles, Cooke+2011b) and for the C- 
enhanced DLA at z abs = 2.34 (square, Cooke+2011a). The color 
shaded areas show the results of the model for the three different 
populations of DLAs predicted to exist at z = 2.34 (see Fig. 3). 

sumably ionized. Moreover, the [Fe/H] value ([C/Fe]) has 
to be interpreted as a lower (upper) limit, since the contri- 
bution by SN type la may be significant in star-rich DLAs 
(Calura, Matteucci & Vladilo 2002). The same interpreta- 
tion has to be applied to [Fe/H] measurements, as Fe is 
depleted onto dust grains (Noterdaeme et al. 2008). Barred 
these limitations, included the unknown number of surviv- 
ing satellites, we note that the apparent concentration of 
[Fe/H]> —2 DLAs in a small region of the Nhi — Fe plane 
(Fig. 4) implies that these absorbers are only partially rep- 
resentative of the progenitors of MW-like systems (Pontzen 
et al. 2008). 



5 DISCUSSION 

The recently discovered z abs ~ 2.34 C-enhanced DLA per- 
tains to a new class of systems, dubbed PopIII DLAs, host- 
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ing the first stars along with second generations of long- 
living stars. These systems are associated to fib-cooling 
minihalos that virialize at z > 8 in neutral and primor- 
dial regions of the MW environment, and only experience 
a short period of SF before the increasing LW background 
burns out their H2. Once "sterilized" these minihalos pas- 
sively evolve as an inert gas reservoir, M g ^ 10 5 ' 5 — 1O 6,3 M0, 
retaining the imprint of the stellar generations they hosted. 
Their gas temperature, T g , is regulated by the balance be- 
tween molecular/metal radiative cooling and photo-heating 
by the external ionizing radiation (Black 1981): 

nni A(T g ,Z) = {e)K ph (2) 

where uhi is the gas density, A(T g ,Z) the gas cooling 
rate (Maio et al. 2007), (e) « 20 eV the mean UV 
background photon energy, and K v h the optically-thick 
H 1 photo-ionization rate (Abel & Mo 1998). K ph « 4.1 x 
10- 15 J2i(AWl0 19 cm- 2 r'V 1 with P = 1.6 and J 21 = 1. 
We find that [Fe/H]> -6 minihalos have T g « (40 - 
1100)K < T V i r , implying that these systems can likely sur- 
vive photo-evaporation thanks to self-shielding and metal 
cooling. Their gas temperature increases with decreasing 
metallicity, in agreement with the results by Kanekar et al. 
(2009). The C-enhanced DLA has T g « 70 K while in 
[Fe/H]w -5 (Z « 1O~ 4 Z ) systems T g « 1000 K. 

The peculiar abundance pattern observed in the z abs « 
2.34 DLA does not result from Z = faint SNe, but 
rather from the enrichment by low-metallicity SNII and 
AGB stars, which may start to form as soon as Z > Z cr = 
10~ 5±1 Z Q . While SNII nucleosynthetic products are mostly 
lost in winds, AGB metals are retained in the ISM, caus- 
ing a dramatic increase of [C/Fe]. The chemical evolution of 
[Fe/H]> —5 DLAs is independent on the assumed yields and 
IMF of PopIII stars, confirming the role of ordinary PopII 
stars in driving the enrichment of VMP systems (SSF07), 
recently emphasized by the detection of a Z < 5 x 10~ 5 Z Q 
star with a "normal" chemical abundance pattern (Caffau 
et al. 2011). If Z cr < 10 -4 Z@, as suggested by the exis- 
tence of such star, the PopIII-PopII transition would be even 
quicker. The mass of relic stars in PopIII DLAs is found 
to be M* ?a 10 2_4 M©, implying that they are the gas-rich 
counterpart of the faintest dSphs. 

As stated, the C-enhanced DLA is a minihalo. However, 
we cannot exclude a different interpretation, in which such 
absorber might be a newly formed halo virializing at z « 2.3 
from a rare, metal-free region of the Inter Galactic Medium, 
and actively forming PopIII stars. Since the total mass of the 
star-forming DLA is M h > 10 9 ' 5 M© « M sf (z = 2.3), strong 
feedback is required to substantially remove the initial gas 
mass and match the observed Nm- By determining the dark 
matter content and SF rate of the C-enhanced DLA, it would 
be possible to disentangle these two pictures. We finally note 
that our simple semi-analytical model, which holds similar- 
ities with that proposed by Abel & Mo (1998) for Lyman 
Limit Systems, prevent us from making specific predictions 
on the number of DLAs at z = 2.34. However, the relative 
contribution of C-enhanced DLAs to the total population is 
expected to be extremely low, ?a 0.01%, thus explaining the 
rarity of these systems. 
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